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Abstract 

A new set of laws of motion for turbulent jets propagating in an in- 
tergalactic medium characterized by a decreasing density is derived by 
applying conservation of momentum flux both in the classical and rela¬ 
tivistic framework. Two characteristic features of radio-galaxies, such as 
oscillations and curvature, are modeled by a classical helicoidal jet. A 
third feature of a radio-galaxy, the appearance of knots, is explained as 
an effect due to the theory of images. 
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1 Introduction 


The study of extra-galactic jets started with the observations of NGC 4486 
(M87) wher e ‘a c urious straight ray lies in a sharp gap in the nebulosity ... 
see ICurtisI ('1 9181 '). More recently the analysis of the radio images of extra- 
galactic jets suggests some interesting problems of physics that should be solved. 
One such problem is that of the decreas e in the velocity a s a fu nc tion o f 


the distance from t he parent nucleu s, see Laing and Bridlel ( 2002allb . 2004 1: 


lLaing et all ( 2006h : Laing and Bridlel (2014). Various mechanisms have been 
suggested for the jet’s deceleration, we report some of them: an interaction 
between a relativi stic jet and t he th ermal radiation from an Active Galactic 
Nuclei (AGN), see Melia and Konig] (Il989h : the incorporation of mass in two- 


dimensional relativistic jets, see [Bowman et al.1 (119961 ): continuous injection of 


plasma a t the base of t he jet and dissipation at some distance from the central 


core, see Wang et al 
bility, 


see 


(200if 


a rotatio n-induced Rayleigh Taylor type insta- 
Meliani and Reopens (120091): and the loa ding of stellar mass pro¬ 


duce^ irpelliptical galaxie s, see 


Perucho et al. ( 20141) . Other authors, such as 


lHardcastle and Sakellioul ( 20041) . quote a terminal velocity of 0.3 c and suggest 
that the jet’s velocity is constant. A second problem is that of the physical 
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mechanism which bends the jets origi n ating; in the the head-tail ra dio-galaxies, 
such as NGC 1265 , see 
rtl979i) . or 1638 + 538 


Owen et al. ( 1978 ), or 115 9 + 583, see iBurns et al'l 


see 


Burns and Owen ( 1980l) . The suggested physical 


mechanisms for bending are: trajectories in the independent blob model, see 
.laffe and Perolal ( 1 97.‘>l) : an adiabatic model in which the ben ding is produced 


by the ram pressure of a central galaxy in unifor m motion, seelQ ’Donoghuc et al 
( 19931 ): the ram pressure in cluster mergers, see lSakelliou and MerrifieldT i 2000l) . 
The extragalactic radio-jets are characterized by a radio-luminosity which is ex¬ 
pressed in watts and also the flux of kinetic energy is expressed in watts. The 
conversion of the flux of kinetic energy in radio-lum inos ity through a tu r bulen t 
cascade has became an active field of research, see Bicknell and Melrosel ( 1982 1. 

We now pose the following questions. 


• Can the physics of turbulent jets, which are observed in the laboratory, 
be applied to extra-galactic jets? 

• Is it possible to generalize the physics of turbulent jets to a medium with 
a varying density? 

• Can we extend to the relativistic regime the physics of turbulent jets? 

• Can we explain the curved trajectories of the extra-galactic jets such as 
NGC 1265 ? 

• Is it possible to build the image of a curved turbulent jet which is emitting 
synchrotron radiation? 


In order to answer these questions, we derive the differential equations which 
model the classical and relativistic momentum conservation for a jet in the 
presence of three types of medium, see Sections [ 2 ] and [3] A model for the 
composition of a decreasing jet velocity along the i-axis with constant velocity 
of the hosting galaxy along the y-axis is derived, see Section U The presence of 
oscillations is modeled by a hclicoidal jet and a bent helicoidal jet, see Section 
[5j Section [ 6 ] updates an algorithm which allows building the radio-image of a 
bent helicoidal jet. 


2 Classical turbulent jets 


2.1 Luminosity conversion 

The tota l power, (). release d in a turbulent cascade of the Kolmogorov type is, 
see 


Pelletier and Zaninetti (11984) 1. 

Q ~ 7 khPSt 


(1) 


where 7 kh = 7 adyr is the growth rate of K-H instabilities, p is the density of 
th e matte r r is the local radius of the jet and st is the velocity of sound, see 


Zaninetti ; 201)7;. The total maximum luminosity, L tl which can be obtained 
for the jet in a given region of radius r and length r is 


L t = 7 T r 2 rQ = ^adpsr^r 2 


( 2 ) 
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The mechanical luminosity of the jet is 


r x 3 2 

= -pv nr 


(3) 


and therefore the efficiency of the conversion, \Ti of the total available energy 
into turbulence is 

_ L t _ 2 1 _ 1 

XT L m Xad 7T M 3 ~ M 3 ’ ^ 

where M is the Mach number. It is clear that in a hot jet characterized by 
high values of M the fraction of the total available energy released firstly in the 
turbulence and after in non-thermal particles is a small fraction of the bulk flow 
energy. The assumption made in the following, in which the bulk flow motion 
is treated independently from t he non-thermal emission, is now justif i ed. A 
similar approach can be found in Bicknell and Melrosel ( 1982 ); Bicknelll ( 1984 ). 


2.2 The parameters for a turbulent jet 

The physics of a turbulent jet can be divided into the simple model and the 
complex model. The simple model is characterized by an opening angle a, and 
th e mat t er’s d ensity p is the same inside and outside the jet, see Section 35 in 
Landau ( 1987 1. The complex model is characterized by the turbulent viscosity, 
vt and, as an example, an opening angle that is a function of the turbulent 
viscosity can be derived starting from Eq. (5.104) in Pone i 


a = 2 arctan ^8 (y/2 — 1^ i't) 


(5) 


In the complex model the matter’s density p is included in vt- In both models, 
the temperature and the pressure are absent and we can speak of cold jets. 


2.3 Momentum conservation 


The co nservation of the momentum flux in a ‘turbulent jet’ as outlined in lLan dau 
(1987) requires the perpendicular section to the motion along the Cartesian x- 
axis, A 

A(r) =7r r 2 (6) 

where r is the radius of the jet. The section A at position Xq is 


,a , 


A(x 0 ) = n(x 0 tan(^)) 2 


(7) 


where a is the opening angle and Xq is the initial position on the x-axis. At 
position x we have 

A{x) = 7r(xtan(|)) 2 . (8) 

The conservation of momentum flux states that 


p(x 0 )uoA(xo) = p(x)v(x) 2 A(x) , (9) 

where v(x) is the velocity at position x and Vo(xo) is the velocity at position xo- 
The selected physical units are pc for length and yr for time; with these 
units, the initial velocity vq is expressed in pcyr -1 , 1 yr = 365.25 days. When 
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the initial velocity is expressed in kms -1 , the multiplicative factor 1.02 x 10 -6 
should be applied in order to have the velocity expressed in pcyr . The tests 
are performed o n a t ypical di stance of 15 kpc relative to the jets in 3C31, see 
Figure 2 in ILaing and Bridle ( 2002bh . 


2.4 Constant density 

In the case of a constant density of the intergalactic medium (IGM) along the 
^-direction, the law of conservation of the momentum flux, as given by Eq. m, 
can be written as the differential equation 

(^z(f)) (x(t)f-v o 2 xo 2 = 0 ■ (10) 

The analytical solution of the previous differential equation can be found by 
imposing x = xg at t=0, 


x(t) = \/2 tv 0 x o + x 0 2 


The asymptotic approximation, see 


Olver et al. ( 2010l) . 


is 


( 11 ) 


x{t) ~ VZy/voXoVt. 


The velocity is 


v(t) 


__VqXQ__ 

y/2tVoX 0 + Xg 2 


and its asymptotic approximation 


v(t) 


V2y/vo^x^ 
2 \ft 


( 12 ) 

(13) 


(14) 


The transit time, t tr . necessary to travel a distance x max 
Eq. (fill) 

, _ Xq T X max 

tfr - . 

2 v 0 x 0 


can be derived from 
(15) 


As a numerical example, inserting cco=100 pc, x max = 1510 3 pc = 15 kpc, 
which is the reference value, and a classical initial velocity of vo =10000 km/s ( 
no = 0.0102pc/yr), we obtain t tr = 1.110 s yr. 


2.5 An hyperbolic profile of the density 

The density is assumed to decrease as 


p = Po(—) , (16) 

x 

where po = 0 is the density at x = Xg- The differential equation that models 
the momentum flux is 

£ (t) (t)) — n 0 2 a; 0 = 0 , (17) 
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and its analytical solution is 


z{t) = 


1 (l2tv 0 xo 2 + 8 a; 0 3 ) 


2/3 


x 0 


The asymptotic approximation is 


^ 1 12 2 / 3(^ 0 2 ) 2 / 3 t 2 / 3 

X[) ~ 4 Xo 

The analytical solution for the velocity is 

xqVq 


v(t) = 2 ■ 


\/12tvoxo 2 + 8a;o 3 
and its asymptotic approximation is 

12 2 / 3 u 0 2 / 3 ^ 


v(t) 


6 \ft 


The transit time can be derived from Eq. © 

—2 £ 0 3/2 + 2 x max 3 / 2 


tf.r — 


3 y/XoVQ 


(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 


and with the same parameters as in Section [2.41 we have t tr = 1.19 10 7 yr. 


2.6 An inverse power law profile of the density 

The density is assumed to decrease as 


P = Po(~) 5 , (23) 

x 

where po is the density at x = Xo- The differential equation that models the 
momentum flux is 

(—) v 2 x 2 — Vo = 0 . (24) 

There is no analytical solution, and we simply express the velocity as a function 
of the position, x, 


v(x) 


XqVq 1 

* vW 


(25) 


see Figure [D 


3 Relativistic turbulent jets 

The conservation of the momentum flux in special relativity, SR, in the presence 
of the velocity v along one direction states that 

O(-) 2 1 y2 + p)A(x) = cost , (26) 

C 1 —IT 

c* 
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Figure 1: Classical velocity as a function of the distance from the nucleus when 
xq =100 pc and vq = 10000km/s: (5 = 0 (full line), <5 = 1 (dashes), <5 = 1.2 
(dot-dash-dot-dash) and <5 = 1.4 (dotted). 


where A(x) is the considered area in the direction perpendicular to the motion, 
c is the speed of light, and w is the enthalpy per unit volume 


w = cp + p 


(27) 


where p is the pressure, see Gourgoulhon ((2006) or formula A30 in De Youngl 
( 2002| ). In accordance with the current models of classical turbulent jets, we 
insert p = 0 and the conservation law for relativistic momentum flux is 


(. Pv 2 1 y 2 )A(x) = cost . (28) 

1 —!T 

er 

Our physical units are pc for length and yr for time, and in these units, the 
speed of light is c = 0.306 pc yr -1 . 


3.1 Constant density in SR 

The conservation of the relativistic momentum flux when the density is constant 
can be written as the differential equation 

(s l(t) ) 'C0) 2 (t™ (f)) a ( J - 

-v 0 2 nx 0 2 (tan ~ “r) =0 • ( 29 ) 

An analytical solution of the previous differential equations at the moment of 
writing does not exist but we can provide a power series solution of the form 

x(t) = dp + ci±t + ci 2 1 2 + + • ■ • , (30) 
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see 


4 are 


Tenenbaum and Pollard ( 19631) : Ince ( 2012 ). The coefficients a n up to order 


a 0 = 


xo 

ai = 


Po c 


1 

3 c 2 /3 0 4 ^ 0 — 3 c 2 {3q 2 xq 

&2 — 

6 

x 0 2 


1 4 c 3 

/3o 7 — 7 c 3 /3 0 5 +3 c 3 /3o 3 

a 3 

6 

xi? 


( 31 ) 

where /?o = ^. 

In order to find a numerical solution of the above differential equation we 
isolate the velocity from eq. HMD 


v{x;x 0 = 


and separate the variables 


XoPqC 2 


\Z—c 2 x 2 I3o 2 + c 2 f3o 2 xo 2 + c 2 x 2 


\/—c 2 x 2 l3o 2 + c 2 /3o 2 xo 2 + c 2 x 2 


dx = dt . 


(32) 


(33) 


Jx o X 0 /3 0 C 2 JU 

The integral on the left side of the previous equation has an analytical solution 
and the following non-linear equation is obtained 

AN 


AD 


= t 


(34) 


where 


AN = /3 0 2 x q 2 In o 2 + lx + \J—x 2 /3 0 2 + po 2 x 0 2 + x 2 ^ 


—fio 2 Xo~ 111 ( —0O 2 + 1 + 1 


-/3 0 2 x 0 2 In (x 0 ) + x\J -x 2 p 0 2 + (3 0 2 x 0 2 + x 2 \j-f. 3 0 2 + 1 


~ x q 2 \l —fdo 2 + 1 


and 


(35) 

(36) 


AD = 2 \J—flo 2 + lx 0 /?o • 

Figure [5] reports an example of the above numerical solution. _ 

Mor e details on the case of constant density can be found in Zaninetti 
(|2009bh . 


3.2 A hyperbolic density profile in SR 

The conservation of the relativistic momentum flux in the presence of an hyper¬ 
bolic density profile as given by Eq. (HI is 

x° x (t) (^x (t)\ 7r (tan (a/2)) 2 (l- ^ dt J — ^ 


—vo 2 n xq 2 (tan (a/2)) I 1 — —h- ) =0 . 


vo 


c 

2 \ - 1 
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Figure 2: Non-linear solution as given by Eq. (l34l) (full line) and series solution 
as given by Eq. (1501) (dashed line) when xq =100 pc and 0o =0.999. 


The analytical solution of the previous differential equation is 

(3ct0 o 3 - 3 ctf3 0 - 2 x 0 ) 2 + 2/3 0 2 x 0 


x(t) = 


The asymptotic approximation is 


2A.-2 


1 ^^„V(A 2 -!) 2 p 0 2 x 0 


:i 5 (A 2 -1) (i-') 2/3 


00-1 

2 ^o 4/3 ^9\/0o 2 c 2 (0o 2 - 1 ) 2 v / F t 

+ 9 0 o c ( 0 o 2 -!) 2 

The analytical solution when the velocity is expressed as 0 = v/c is 

QU ^ _ -0o (3 0o 3 ct - 3 0o ct - 2 x 0 ) s/2 

^ * “ 77 7 ^ 2\275 ’ 

(J3 0 o 3 cf - 3 0 o ct - 2 x 0 ) J 

and the asymptotic approximation is 

-0o 2 c (0 O 2 - 1) 


m = 


3 (0o 2 c 2 (0o 2 - l) 2 ) 


2/3 


(38) 


(39) 


(40) 


(41) 


3.3 Inverse power law profile of density in SR 

The conservation of the relativistic momentum flux in the presence of a density 
profile of the inverse power law type as given by Eq. (1231) is 

-xo 5 x~ s+2 /3q' 2 i3' 2 + 0 O 2 0 2 x 0 2 + x 0 s x~ s+2 /3 2 - 0 O 2 x 0 2 

(02 - 1) (0 O 2 - 1) 
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= 0 . (42) 














Figure 3: Relativistic variable j3 as a function of the distance from the nucleus 
when xq =100 pc and /?o =0.9: (5 = 0 (full line), 8=1 (dashes), 8 = 1.2 
(dot-dash-dot-dash) and 8 = 1.4 (dotted). 


This differential equation does not have an analytical solution and the expression 
for /? as a function of the distance is 



The behavior of /3 as a function of the distance for different values of 8 can be 
seen in Figure [3] 

The transit time can be derived from Eq. m 


ttr — 


2 x 0 2 - 2 \Jx 0 (/3 q 2 x 0 - /?o 2 


3/3qc (/3 0 2 - l) x 0 


(44) 


As an example, inserting io=100 pc, x max = 1510 3 pc, and /?o = 0.9, we have 
t tr = 2 10 5 yr. 


4 Classical curved trajectories 


In the presence of a 3D trajectory, r(t), the acceleration is given by 


a = 


ck’ v 2 T 

—T + - N 

at p 


(45) 


where v = = |^f is the magnitude of the velocity, s is the arc-length of 

the trajectory, T and N are the tangential and normal versors to the trajectory, 
and R is the radius of curvature given by 


R = 


,d x 


,d 2 u.n ,drz. 9 
y ds 2> y ds 2> 


- 1/2 


(46) 
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see exercise 5.30 in Spiegel ( 197lh . A second definition of the radius of curvature 
when we have a 2D curve parametrized as y(x) is 


R = 


1 + 


( dy\ 

\ dx J 


d 2 y 
dx 2 


(47) 


see Eq. (12.5) in Granville ( 19111 ). In the presence of a trajectory given by a dis¬ 
crete set of points, the circle of curvature which has radius R c , and coordinates 
of its centre ( x c ,y c ), can be drawn w hen three points are given, see Appendix 


1X1 According to theorem 14.1.1 in Granvillel ( 1911 b the radius of the circle of 


curvature and the radius of curvature are equal. 


4.1 The astronomical data 

As an application, we analyse the first part of the right side of NGC 1265 where 
the distances were evaluated adopting Hq = 50kms~ 1 Mpc^ 1 for the Hubble 
constant and z = 0.0183 for the redshift see Xu et al. (Il999il . At the moment 
of writing, there is a more precise eval uatio n of the Hubble constant, which is 


H q = (69.6 ± 0.7)kms 1 Mpc 1 , see Bennett et al. (120141), which , coupled with 


c = 299792.458-^2 for the velocity of light, see Mohr et all (2 012 1. means 

1" = 382.15pc . (48) 

The trajectory of NGC 1265 was derived as follows 

1. The data were digitized in " from Figure 3 in Owen et al. (1197811 using 


WebPlotDigitizer, a Web based tool to extract data from plots. 
2. The conversion from " to pc was done using Eq. 


Figure |H shows the digitalized 2D trajectory and 8 circles of curvature computed 
according to Eqs (1A.1I IA.2IA.3ll . 

The average radius of the circle of curvature, R Cl for the first 20 kpc of the 
right side of NGC 1265 is 


R c = (9.22 ± 1.03)fcpc . 

4.2 Curved trajectory with constant density 


(49) 


A ballistic trajectory is determined on Earth by gravity and aerodynamic drag. 
Here we assume a non-ballistic trajectory for the radio-galaxy as being due to 
the composition of the jet’s motion along the i-axis with a motion at constant 
velocity, v g , of the parent galaxy along the y-axis. The peculi ar v elocity for 
a gala xy with redshift z can be evaluated using formula (3) in Freeland et al.1 
(20081). Other authors replace the velocity of the p arent galaxy w i th the velocity 
of external winds, see lHardee and Hughesl ( 2003h : LPerkins et al.l ( 2004h . 

When the galaxy’s velocity is expressed in kms -1 , the multiplicative factor 
1.02 x 10“ 6 should be applied in order to have the galaxy’s velocity expressed 
in pcyr^ 1 . The two equations of motion are 


x{t) = \Z2tv 0 x 0 + Xq 2 
y(t ) = yo + v g t , 


(50a) 

(50b) 
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Figure 4: The right side of NGC 1265 , the line with a large width, represents 
the real data as extracted by the author from Figure 3 in Owen, Burns, & 
Rudnick (1978) and 8 circles of curvature. 


where yo is the y-position at t = 0. In order to evaluate the radius of curvature 
in the case of constant density, we extract t from Eq. (150all and we insert it into 
Eq. (I50bl) 


y{x) = 


1 v g ( x 2 - x 0 2 ) 

2 v 0 x 0 


+ yo 


(51) 


The previous formula allows visualizing the trajectory which represents the best 
fit for the right side of NGC 1265 in the case of constant density, see Figure [5] 
The parameters vo and v g are found by minimizing the value of \ 2 defined 
as 

n 

X 2 = ^2(Vi - Vi,th ) 2 , (52) 

i—1 


where yi is the y-value of the itli point in the digitized trajectory and yi^h is 
the theoretical point obtained by inserting the x value of the *th point in Eq. 

m- 

The radius of curvature in the case of constant density as given by Eq. (HU) 


is 


R = 


(■X 2 Vg 2 + V 0 2 X 0 2 ) 3/2 


V 0 2 X0 2 Vg 


(53) 


as a consequence the centripetal acceleration is 



^2 tv 0 x 0 + x 0 2 + {tv g + y 0 ) 2 ) v 0 2 x 0 2 
(v 0 2 x 0 2 + v g 2 (2 tv 0 x 0 + ir 0 2 )) 3/2 


-N 


(54) 
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Figure 5: Real trajectory of NGC 1265 (big stars) and theoretical curved tra¬ 
jectory (full line) in the case of constant density. The parameters are Xo=100pc, 
?/o=0 pc, vq =9620 krn/s and v g = 124 km/s. 


4.3 Curved trajectory with hyperbolic density 


In the case of an hyperbolic density profile we can isolate the time in Eq. m 
and insert it into Eq. (I50bl) . The resulting trajectory, which is independent of 
the time, is 


y{x) 


2 v g (x 0 3/2 - x 3 / 2 ) 

3 y/XQVQ 


(55) 


The trajectory which represents the best fit for the right side of NGC 1265 in 
the case of an hyperbolic density profile is shown in Figure [6] 

The radius of curvature in the case of an hyperbolic profile of density as 
given by Eq. HMD is 


/ 2 i 2 \3/2 I — 

R _ 2 (xV g ~ + V Q-ttoj y/S 

V 0 2 X 0 Vg 

and the centripetal acceleration is 


(A + (12 tvoXo 2 + 8x 0 3 ) 4/3 + 16 xo 2 yo 2 ^ v 0 2 x 0 '/iv g 




-N 


32 x 0 2 ( V 0 2 X 0 + 1 /4 ^ 2 (12tW+8x 0 3) 2 /3 j 3 / 2 ^ (12tv 0 x 0 2+8x^ X 


3 - 12/3 


where 

A = 16t 2 v g 2 xo 2 + 32tv g xo 2 yo 


(56) 


(57) 


(58) 
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Figure 6: Real trajectory of NGC 1265 (big stars) and theoretical curved tra¬ 
jectory (full line) in the case of an hyperbolic profile of density. The parameters 
are £o=100 pc, yo=0 pc, Vq =3660 km/s and v g =54.74 km/s. 

5 The helicoidal trajectory 

The circular helix is 


x = a cos(t) (59a) 

y = asin(t) (59b) 

z = bt (59c) 

where a is the radius and the pitch is 2i r 6, see Lipschut3 ( 19691 ). The radius of 
curvature of the circular helix is 


b 2 


R = 


The arc-length of the helix, s, as a function of time is 


s = sj a 2 + b 2 t . 

The helicoidal jet in the case of constant density is 


x = \/2tv 0 Xo -(- xo 2 
y = at cos (f l p t) 
z = at sin (f l p t) , 


(60) 


(61) 


(62a) 

(62b) 

(62c) 


where the x—t relationship is given by Eq. cn fl p is the angular velocity, and 
the radius of the helix grows linearly with time according to a linear relationship 
given by the parameter a. The pitch of the helicoidal jet is 


. tt v 0 x 0 
4 — -1- x 0 2 


rip 

The radius of curvature of the helicoidal jet is 

NH 

R = 


DH 


(63) 


(64) 
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Figure 7: The radius of curvature of the helicoidal jet versus time when the 
parameters are a;o=100pc, yo=Opc, and Vq =9620 km/s. 


where 

NH = (2 a 2 t 3 fl p 2 vo + a 2 t 2 tl p 2 xo + 2 a 2 tv o + a 2 x o + vq 2 Xo) 3 ^ 2 , (65) 


DH = 

a ^8 a 2 t 7 n p 6 v o 3 + 12 a 2 t 6 fl p 6 v o 2 xq + 6 a 2 t 5 fl p 6 voXo 2 + a 2 t 4 fl p 6 xo 3 
+32 a 2 t 5 Q p 4 vo 3 + 48 a 2 t 4 fl p 4 Vq 2 Xo + 4t 4 fi p 4 vo 4 Xo + 24 a 2 t 3 fl p 4 voXo 2 
+4t 3 £i p 4 vo 3 Xo 2 + 4a 2 t 2 fl p 4 xo 3 + t 2 tt p 4 v o 2 Xo 3 + 32 a 2 t 3 Q p 2 vo 3 

+48 a 2 t 2 fl p 2 vo 2 xo + 21 t 2 fl p 2 vo 4 xo + 24 a 2 tfl p 2 voxo 2 

0/2 

+18tf2 p 2 uo 3 a;o 2 + 4a 2 f2 p 2 £o 3 + 4 VL p 2 vq 2 xq 3 + vq 4 xq ) . (66) 

Figure [7] shows the radius of curvature of the helicoidal jet in kpc as a function 
of time in units of 10' years. 

The bent helicoidal jet in the case of constant density is 

x = y/ 2 tvoxo + xq 2 (67a) 

y = at cos (flpt) (67b) 

z = at sin (fl p t) + zq + v g t . (67c) 


A 3D display of the bent helicoidal jet is shown in Figure 0 where the choice 
of the Euler angles, which define the observer, corresponds to the astronomical 
observations. 

Another choice of the Euler angles produces a different projected surface, 
see Figure [9] 

The radius of curvature of the bent helicoidal jet has a complicated expres¬ 
sion and we only present its numerical behavior, see Figure flTil 

The arc-length of the bent helicoidal jet is 

,= l 'M* t ■ m 
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Figure 8: Continuous 3D surface of NGC 1265 : the three Eulerian angles 
characterizing the point of view are <f>= 0 °, 0= 90 °, and 5'= 0 °. The 
physical parameters are xo=100 pc, yo=0 pc, vq =9620 km/s, v g =124 km/s, 
fl p = 6 10 -8 7r r e ^. s and alpha = 5°. The axes are expressed in kpc units. 



Figure 9: Continuous 3D surface of NGC 1265 : the three Eulerian angles 
characterizing the point of view are <£= 70 °, 0= 20 °, and 4^= 10 °. The 
physical parameters are the same as in Figure [8] and the axes are expressed in 
kpc units. 
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Figure 10: The radius of curvature of the bent helicoidal jet versus time. The 
physical parameters are the same as Figure [5] 


where 

NS = 2 a 2 flp 2 t 3 v o + a 2 fl p 2 t 2 xo + 4 cos (A p t) aCl p t 2 vov g + 2 a 2 tv o 
+2 cos (Q p t) aflptVgXo + a 2 x o + 4 sin (f l p t) atvov g 
+2 sin (fl p t) avgXo + 2 tvov g 2 + vo 2 xo + v g 2 xo , (69) 

and 

DS = 2tvo + xo . (70) 

The previous integral does not have an analytical expression and the numerical 
integration between 0 and 10 8 years, the other parameters as in Figure [10] gives 
s = 22.2 pc. 


6 The theory of images 

In the case of an optically thin layer, the emissivity is proportional to the number 
density, C, 

ju C = KC{s) , (71) 

where is the emission coefficient and K is a constant. This can be the case 
for synchrotron radiation in the presence of an isotropic distribution of electrons 
with a power law distribution in energy, N(E), 


N(E)dE = K S E~ 


(72) 


where K s and 7 are constants, see lZaninetti ( 2010ll for more details. In the case 
of constant number density 

jv(, = K C s , (73) 
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Figure 11: The source is represented through a circular section perpendicular 
to the jet axis. The observer is situated along the y-direction, one line of sight 
is indicated and the angle /3 is clearly indicated. 


where s is the length of the relevant line of sight, this formula is extremely simple 
and allows building the image in the appropriate geometrical environment. Two 
analytical results outline the theoretical framework that should be verified by 
the simulation. The first one analyses the behavior of the intensity or brightness 
along a jet when the distance from the origin, x, is fixed. We assume that the 
number density C is constant in a cross section of radius a and then falls to 0, 
see Figure mi 

The length of sight, when the observer is situated at the infinity of the y-axis, 
is the locus parallel to the y-axis which crosses the position x in a Cartesian 
x — y plane and terminates at the external circle of radius a. The locus length 
is 


l a b = 2 x ( sja 2 — x 2 ) ; 0 < x < a . (74) 

When the number density C is constant on a cylinder of radius a, the intensity 
or brightness of radiation is 

loa = C x 2 x (\/a 2 — x 2 ) ; 0 < x < a , (75) 

or 

loa = C x 2 x a x sin(/3) ; ~ < /3 < ^ , (76) 

which can be called the ‘trigonometrical law’ for the intensity or brightness. 
The second analytical result can be derived when the curved shape of a jet of 
finite cross section is parametrized by a toroidal helix. The toroidal helix has 
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Figure 12: The section of one-fourth of a toroidal helix is represented through 
a circle of radius R — r and a bigger circle of radius R + r. The observer 
is situated along the y-direction at infinity and the line of sight of maximum 
length is indicated. 


the following parametric equations: 


x = cos(a) • (R + r ■ cos(0)) 

y = sin (a) • (R + r ■ cos(9)) (77) 

Tl T 

z = r • sin(0) + a-— , 

Z7T 

(78) 


where 0 £ [0, 27t), a £ [0, 2n), nr is the distance along the 2 -axis after an angle 
a = 27r, and n is an integer. We now analyse the case in which 

a < — , (79) 

n 

where the right-hand side is the value of the angle after which the toroidal 
helix has advanced by r along the 2 -direction. Figure [T3] shows a section in the 
middle of the toroidal helix 2 = 0, from which is possible to see that the dotted 
line presents the longest line of sight, l max , which starts from x = R when the 
observer is at the infinity of the y-axis. 

The shortest line of sight is 2r. The maximum enhancement in the presence 
of a constant number density, e, is 


max 



A simple geometrical demonstration gives 


e 


f V2 R + r 
2 Vr 


(80) 


(81) 
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The relationship between the radius of curvature, R , and the radius of the 
toroidal helix, r, which produces an enhancement e in the intensity or brightness 
is 



(82) 


We now outline how it is possible to build a radio-image. The number 
density is stored on a 3D grid k) where i,j and k are indices varying 

from 1 to pixels. The orientation of the object is characteriz ed by t he Euler 
angles (<b, 0,4') and therefore by a 3 x 3 rotation matrix, E, see Goldstein et al.1 
( 2002h . The grid k ) is then rotated according to the chosen Euler angles. 

The intensity map is obtained by summing the points of the rotated images along 
a particular direction, and the intensity is 


= ^Asx M(i,j,k) (83) 

k 

optically thin layer linear case , 


where As is the spatial interval between the various values of intensity and the 
sum is performed over the interval of existence of the index k. In this grid the 
little squares that are characterized by the position of the indexes i.j correspond 
to a different line of sight. 

The effect of the insertion of a threshold intensity, It r , given by the observa¬ 
tional techniques, is now analysed. The threshold intensity can be parametrized 
to Imaxi the maximum value of intensity characterizing the map, 


Itr — 


Imax 

factor 


(84) 


where the parameter factor is greater than one. 

A map for the emissivity of a jet can be built by employing the following 
algorithm: 


• A great number of points, for example, one million, is inserted into a bent 
helicoidal jet with given characteristic parameters. 

• This set of points is rotated according to the three Euler angles which 
identify the observer’s point of view. 

• The 2D matrix which represents the intensity of the non-thermal emission 
is built according to procedure (1M|) . 


Figure [13] shows the projected points. The two analytical results, the enhance¬ 
ment of the intensity in the central line and the knot structure due to the 
curvature along the line of sight, are clearly visible in Figure 1141 which shows 
the intensity of a non-thermal map of NGC 1265 . 

A first comparison of the previous figure can be done with Figure 1 in 


Alo v et ah (120031) where a 3D relativistic hydrodynamic simulation for the 
precessing beam was carried out. A second comparison can be done with 
the ps e udo-synchrotr o n inte nsities visible in the figures of lHardee and Hughesl 
(2003j); Hardee et ajJ ([2005) where the knots in the radio-jets were simulated 
in the framework of helical relativistic instabilities. A third comparison can be 
done with Figure 3 in 'Laing and Bridle! ( 20141) where the observed and model 
total-intensity images are reported for 15 radio-galaxies without bending effects. 
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Figure 13: Continuous 3D trajectory of NGC 1265 : the three Eulerian angles 
characterizing the point of view are $= 0 °, 0= 90 °, and \E'= 0 °. The physical 
parameters are the same as Figure [8] 


minimum= 0.0139 maximum= 1.0000 
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Figure 14: Theoretical 2D map of the surface brightness of the emission of 
NGC 1265 with basic parameters as in Figure fl3l The integration is performed 
on a cubic grid of 1000 3 pixels and factor = 4. 
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7 Conclusions 


Classical turbulence: Turbulent jets are usually modeled by a temporal evo¬ 
lution with density equal to that of the surrounding medium. Here, in order 
to cover the astrophysical applications, we considered an hyperbolic profile of 
density, see the solution ED and an asymptotic solution flUD . The case of a 
density which follows an inverse power law is limited to the derivation of the 
velocity, see Eq. HO). This inverse power law case allows matching a wide 
variety of astrophysical situations, as an example, when 6 = 2 , the velocity does 
not decreases with distance. 

This is obvious from the conservation equation, pvA = constant. If A goes 
as r 2 and p goes with the inverse of the square of distance in a conical jet, v is 
constant. 

Relativistic turbulence: The conservation of the relativistic momentum 
flux for turbulent jets is here analysed in three cases. The first case is that with 
a surrounding medium having constant density, where the analytical result is 
limited to a series expansion for the solution, see Eq. (ED. The second case is 
that of an hyperbolic density decrease for the surrounding medium, for which 
we derived an analytical solution see Eq. (1551) and an asymptotic solution, see 
Eq. ED- The third case is that where the surrounding density decreases with 
a power law behavior: the analytical result is limited to the velocity-distance 
relationship, see Eq. ED- 

The curvature of the jet: The composition of the velocity is discussed in 
the light of the radius of curvature and the standard mathematical definition is 
given, see Eq. ED- The astronomical counterpart is the radius of the circle of 
curvature which is shown in Figure [I] for NGC1265 . Two analytical solutions 
are given for the curved trajectory, see Eqs m and (l50bl> . in the case of 
constant density and Eq. (l55l) for an hyperbolic decrease in density. A careful 
analysis of Figures [5] and [6] does not show an accurate coincidence between the 
predictions of the model and the digitized trajectory. In this case, its y 2 can be 
lowered by introducing other effects such as photon losses due to synchrotron 
radiation or another density profile. 

Helicoidal bent jet: The precession is here modeled by an helicoidal jet, 
see Eq. [62j which has a pitch and radius of curvature given by Eqs (1M1) and 
ED- The composition with the velocity of the host galaxy allows modelling 
the a helicoidal jet in which the radius of curvature can be visualized only 
numerically, see Figure El The arc-length has a complicated expression, which 
is given by Eq. (USD- 

Theoretical radio maps: The radio image of an extragalactic radio source 
is built adopting: (i) a uniform number density of synchrotron emitters over the 
entire jet, (ii) the thin layer approximation. The two theoretical effects of central 
brightening, see Eq. ED- and intensity brightening due to the curvature of the 
emitting region, see Eq. ED- are both visible in the simulation of NGC 1265 , 
see Figure [E 

New analytical results can be obtained expressing the helicoidal trajectory 
in cylindrical coordinates or analysing the ‘valley on the top’ effect which is due 
to the velocity profile in the direction perpendicular to the motion, see Zaninetti 
(2009a). The previously cited papers solve some of the problems connected with 
the radio-jets but leave other problems open, and Table [T| reports their status. 


21 








Table 1: Synoptic table of the assumptions of two papers and here 


Problem 

LainghBridle 2002 

Hardee et al.2005 

this paper 

distance along the jet 

numerical 

not 

analytical 

velocity along the jet 

numerical 

not 

analytical 

knots 

not 

instabilities 

image theory 


A The circle of curvature 

Once three points are selected on a curve, (x\,yi), (x 2 , 2/2)1 (£3,2/3), the circle of 
curvature has radius 

■ <-> 

where 

NR = | {x 2 2 - 2 x 2 x 3 + x 3 + y 2 2 -2 y 2 y 3 + y 3 2 ) x 

(xi 2 - 2 X\X 3 + x 3 2 + y 2 - 22 / 12/3 + y 3 2 ) x 

(1/2) 

(x! 2 -2xix 2 +x 2 2 + yi 2 -2y 1 y 2 + y 2 2 ) l , 

DR = 2 xiy 2 - 2 xiy 3 - 2 x 2 y 3 + 2 x 2 y 3 + 2 x 3 yi - 2 x 3 y 2 . 

The two coordinates of the centre of the circle, (x c ,y c ), are 

NXC 


X c = 


DXC 


where 


and 


Vc = 


NYC 

DYC 


where 


NYC = -x 2 x 2 + x\ 2 x 3 + xix 2 2 - x\x 3 2 + xiy 2 2 - xry 3 2 
-X 2 x 3 + £ 2 £ 3 2 - x 2 y\ 2 + x 2 y 3 2 + x 3 y 3 2 - x 3 y 2 , 

DY 0 = 2 xiy 2 - 2 x 3 y 3 - 2 x 2 yi + 2 x 2 y 3 + 2 x 3 y 3 - 2 x 3 y 2 . 


(A.2) 


NXC = xi 2 y 2 - x\ 2 y 3 - x 2 2 y! + x 2 2 y 3 + x 3 2 yi - x 3 y 2 + y 2 y 2 
-yi 2 V 3 - 2/12/2 2 + 2/12/3 2 + y2 2 y 3 - 2/22/3 2 , 

DXC = 2 xiy 2 - 2 x 3 y 3 - 2 x 2 yi + 2 x 2 y 3 + 2 x 3 y 3 - 2 x 3 y 2 , 


(A.3) 
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